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Edge effects in a frustrated Josephson junction array with modulated couplings 
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A square array of Josephson junctions with modulated strength in a magnetic field with half 
a flux quantum per plaquette is studied by analytic arguments and dynamical simulations. The 
modulation is such that alternate columns of junctions are of different strength to the rest. Previous 
work has shown that this system undergoes an XY followed by an Ising-like vortex lattice disordering 
transition at a lower temperature. We argue that resistance measurements are a possible probe of the 
vortex lattice disordering transition as the linear resistance Rl(T) ~ A(T)/L with A(T) oc (T — T c i) 
at intermediate temperatures T c xy > T > T c i due to dissipation at the array edges for a particular 
geometry and vanishes for other geometries. Extensive dynamical simulations are performed which 
support the qualitative physical arguments. 
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There has been a lot of interest in arrays of supercon- 
ducting grains coupled by Josephson junctions. In the 
absence of a magnetic field perpendicular to the plane 
of the array, such a system is rather well described by 
a classical two dimensional rXY model and experimen- 
tal transport measuixmcntsEI are well described by the 
dynamical extension! 2 ! of the static theoryHa. The most 
recent technique employed to measure truly equilibrium 
properties of a superconducting array is a magnetic flux 
noise measurementEj which does not require an imposed 
current. Earlier experimental studies measured the yolU 
age V due to an applied zero frequency current iHiLjj 
and the frequency dependent impedance Z(lo) by two coil 
mutual inductance techniques!!!. The agreement between 
theory and experiment is quite good. In the presence 
of a magnetic field normal to the array, the situation is 
not so clear as, even in the simplest case of half a flux 
quantum per elementary plaquette / = 1/2 on a square 
lattice, the system is much more complicated and less 
well understood. Even in the absence of disorder, which 
is inevitably present in an experimental system!! 2 !, there 
are two types of competing order when / = 1/2: a dis- 
crete Z2 symmetry of the ground state of the vortex lat- 
tice and the U(l) symmetry of the superconducting order 
parameter. In an isotropic system where all the junction 
strengths are the same the phase coherence of the su- 
perconducting order parameter must be destroyed at a 
lower temperature or, at best, the same,temperature as 
the discrete order of the-s^ortex latticeO. The most re- 
cent simulations on thiscJ'Ej agree that the phase order 
is destroyed by a KT transition at a very slightly lower 
temperature than the discrete order which undergoes a 
transition in the Ising universality class. 

An interesting variant of the isotropic souare array 
with / = 1/2 was proposed by Berge et aO in which 
alternate columns of junctions are of different strength 
to the rest. This model is described by the Hamiltonian 
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where K v v i = K = J/ksT on all nearest neighbor bonds 



except on alternating columns where K r 
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A(r, r') = ^ J r A • dl where A is the vector poten- 
tial of the external magnetic field and $0 is the quan- 
tum of flux. The frustration of plaquette R is /(R) = 
Sdr ^( r ' r ')/2 7r = 1/2 where the sum is over the bonds 
< r, r' > surrounding the plaquette. It is a gauge in- 
variant definition and a convenient gauge is the Landau 
gauge in which A(r, r') = n on every second column of 
bonds and A(r, r') = otherwise. Such an array is exper- 
imentally realizable by varying the area of the appropri- 
ate junctions to obtain alternating columps of strength 
r\K and is currently being investigated!!-! The. mndrl. 
described by Eq.(||) has been studied numericallyHu3'E£l 
and it was found that, by varying the anisotropy parame- 
ter 77, the order of the XY and Ising transition is reversed. 
In the isotropic system at 77 = 1 the XY transition oc- 
curs at a lower temperature than the Ising transition but 
when \r\— 1| is sufficiently large the Ising is at a lower tem- 
perature than the XY transition. The early simulations 
were unable..ta separate the transitions for \rj — 1| w but 
later world 1 Tr 5 ! provided convincing evidence for two sep- 
arate transitions with T c j > T c xy- Another generalized 
version of thp_square array at / = 1/2 has been consid- 
ered recentlyE] in which the couplings K r>r +f in Eq. (|l|) 
are modulated in the x and y directions simultaneously, 
in a zig-zag pattern. In this case, varying r\ in the range 
Tj > 0.6 does not lead to a change in the order of the Z2 
and U(l) transitions and they critical properties remain 
the same as for the isotropic square array. 

In a two dimensional plane of superconductor in a mag- 
netic field normal to the plane a vortex lattice is formed. 
If a vortex is subjected to a force due to an external 
applied current, it will move in a direction perpendicu- 
lar to the current and this creates a voltage V ~ nfl 
where nj is the density of free vortices, which implies a 
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finite linear resistance. This is regarded as signalling the 
destruction of superconductivity. An array of Joseph- 
son junctions with no disorder in a magnetic field with 
a frustration / = p/q has a ground state which is a vor- 
tex lattice commensurate with the underlying lattice and 
is pinned and superconducting. As the temperature is 
raised, the vortex lattice will melt or become a floating 
phase which is not pinned. In cither case, one would 
naively expect that the system becomes non supercon- 
ducting as the vortices will move and induce a voltage 
when an external current is applied. In particular, a 
square, array with / = 1/2 and r\ > 1/3 has a ground 
stateEj which is a vortex lattice with a vortex in every 
second plaquette commensurate with the underlying lat- 
tice. Since the vortex lattice is pinned by commensurabil- 
ity effects, the system will be superconducting and should 
become non-superconducting when the lattice melts. In 
the anisotropic array of Eq. (|l|) , the vortex lattice melts by 
an Ising transition at T = T c j but XY ordeaj-feupercon- 
ductivity) persists to higher temperaturescj'Eia although 
the vortex lattice has melted, which seems to contradict 
the standard folklore that, when the vortex lattice melts, 
superconductivity disappears. It seems to us that there 
are three possibilities to reconcile the equilibrium behav- 
ior of the anisotropic system of Eq.(|l|) with these qual- 
itative arguments about the transport properties with 
a small applied current I at intermediate temperature 
T c j < T < T c xy where the vortex lattice is melted: (1) 
the present understanding of the effect of the melting of 
a vortex lattice on superconductivity is wrong; (2) the 
equilibrium calculations at / = have no relation to the 
dynamics as / — > 0; (3) the properties of the melted 
vortex lattice in the anisotropic array are special and the 
qualitative argument that superconductivity is destroyed 
by the melting does not apply. 

In this paper, we take the view that the system of 
Eq.(|l]) is special and that scenario (3) is the explana- 
tion and also address the question of the signature (if 
any) of the low temperature Ising transition in experi- 
mental measurementstH on an anisotropic array with col- 
umn modulation. Two coil mutual inductance experi- 
ments measure the dynamical impedance Z(ui) which is 
proportional to the inverse of the helicity modulusEJ 7. 
At the Ising .transition, this has a harmless singularity of 
the fornxSEj 7 ~ t In t where t = T/T c j — 1 which implies 
that the impedance will also not show any divergence as 
observable signature. The implications for the flux noise 
spectrum have not been worked out and it would be of 
some interest to do this. There is one possibility which is 
the linear resistance which we discuss in the rest of this 
paper. We study by qualitative analytical methods and 
by numerical simulations the zero frequency IV charac- 
teristics of an anisotropic array and show that the onset 
of linear resistance occurs at the low temperature Ising 
transition, which may be detectable by experiment. This 
is an edge effect and, in an L x L array, the linear resis- 
tance R L {T) ~ A{T)/L with A{T) > when T > T cI . 
We begin by defining the model and argue that the onset 
of linear resistance is due to the formation and growth of 



domains of reversed chiral order at the edges of the array. 
The geometry of the array and the direction of the ap- 
plied current is important, and we first define the various 
geometries and boundary conditions and what is to be ex- 
pected on physical grounds in each case. We then discuss 
the results of numerical simulations and show that these 
are consistent with our analytical arguments. 

The low temperature Ising transition may be described 
in terms of the proliferation of domain walls separating 
regions of different ground states. A ground state config- 
ur ation of Eq.(g) may be represented in the Coulomb gas 
representation by a set of unit charges in a-checkerboard 
pattern on half the sites of the dual latticetLa. The impor- 
tant low energy excitations which destroy the Ising order 
are domain walls between regions of different Ising or- 
der and these domain walls lie along the bonds of the 
original lattice. It is well known that there are frac- 
tionaixharges q — ±1/4 at the corners of these domain 
wallstJcj. One may argue that such domains with a 
net integer charge formed of corner charges of the same 
sign are the excitations which undergo a KT unbinding 
transition at the same temperature as-.t hq . Is ing transi- 
tion. However, numerical simulationsEJ^tSlla show that 
the XY and Ising transitions are separated for almost all 
values of the anisotropy parameter rj which implies that 
the fractional corner charges and the integer charges do 
not screen each other. In the case of interest here, \r) — 1| 
sufficiently large, one can understand this by consider- 
ing the domaixL. walls at T c j < T < T c xy- As shown by 
Eikmans et al£3, the domain walls between regions of dif- 
ferent Ising ground states lie along the strong bonds since 
the energy of these is less than those on the weak bonds. 
In fact, they find numerically that there is a negligible 
density of domain walls on weak bonds for T < T c xy 
but the density on the strong bonds grows rapidly when 
T > T c i. Thus, these domains carry zero net charge 
as the fractional corner charges must alternate between 
q = ±1/4. For a domain to have an integer charge, one 
of the vertical walls must lie along a weak bond which 
costs a lot of free energy. The unit cell describing the 
Ising degree of freedom may be taken to be four adja- 
cent elementary plaquettes bounded by strong bonds as 
shown in Fig. (1) and at T > T c j domains with an inte- 
ger number of these unit cells will proliferate. Since these 
domains carry zero net charge and zero dipole moment, 
they cannot contribute to the linear resistance which is 
proportional to n/ , the concentration of thermally ex- 
cited free charges (vortices)H. Because of the zero dipole 
moment, the mechanism of stretching of domains-by the 
applied current I suggested by Mon and TeiteEIl as a 
contribution to the non linear resistance for T < T c j of 
the form V ~ Iexp(-f d (T)/IkT), where f d (T) is the 
domain wall free energy, is not a dominant effect for the 
anisotropic arrays considered in this paper. 
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FIG. 1. (a) States of opposite Ising order of 2 x 2 unit cells. 
Solid (dotted) lines denote strong(weak) bonds. Unit charges 
are represented by x's. (b) A domain of opposite Ising order. 
The domain wall is the dash-dotted line. Filled(unfilled) dots 
at domain wall corners are ±1/4 fractional charges. 



the strong bonds are at the edges as in Fig. (2b) as 
there are no net fractional charges at the ends of these 
edge domains and these will behave just like the domains 
in the bulk of Fig. (lb) which carry no net charge. The 
situation for arrays with weak bonds at the vertical edges 
is very different and this is shown in Fig. (2a). In this 
geometry, one can regard these edges as splitting the unit 
cells of Fig. (2a) in half thus forcing fractional charges of 
opposite signs at the two ends of the domains. Since T > 
T c j, the domain walls may be regarded as having melted 
and having zero line tension as the edge undergoes an 
ordinary transition slaved to the bulk transitionEl Along 
the edges of the array of linear size L, these charges are 
effectively free unbound charges of concentration n/ ~ 
1/L£||. This implies that in an array with weak bonds at 
the vertical edges as in Fig. (2a) there will be a linear 
resistance 

R L {T)~l/L^~L-\T-T cI ) (3) 

provided there are free boundary conditions at these 
edges. This can be realized by uniform current injection. 
Other methods of current injection such as injection from 
superconducting busbars will suppress this effect as the 
busbars repel the vortices (charges) from the ed ge s a nd 
the linear resistance will be reduced to 0(L~ nK )E3o as 
in arrays with strong bonds at the vertical edges. 



We consider a 2N X x 2N y lattice of elementary plaque- 
ttes corresponding to a (2N X + 1) x (2N y + 1) lattice of 
superconducting grains coupled by weak links. There are 
two geometries to consider (a) columns of weak bonds at 
the two vertical edges and (b) columns of strong bonds 
at the edges as sketched in Figs. (2a) and (2b). In case 
(b), there are N x x N y complete unit cells and the current 
is injected uniformly in the x direction perpendicular to 
the weak bonds. A linear resistance is due to the cur- 
rent driving the pre-existing thermally excited charges 
across the system but such charges do not exist in the 
bulk because, near the Ising transition, only neutral do- 
mains are present as the domain walls are constrained 
to lie along the strong bondstS. One may expect that 
free, unbound fractional charges exist near the edges of 
the array because of domain formation at the edges. We 
expect that for T > T c ] the system will consist of a set 
of domains of opposite Ising order of size determined by 
the bulk correlation length ^ ~ A\>t~ Vh in the bulk and 
at the edges for free boundary conditions by the_anrface 
correlation length £|| ~ A\\t~ v w with v h = u\\ = ltl The 
amplitudes Af, and A|| are nonuniversal but their ratio is 
universal and is givetuin terms of critical exponents from 
conformal invarianceEH by 

A|,/^ 6 = 776/T7H = 1/4 (2) 

with the bulk and edge exponents@~0 rjb = 1/4 and 
?7|| = 1 which implies that the density of edge domains 
is larger than that in the bulk of the array. This is of 
no significance for the linear resistance in arrays where 
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FIG. 2. (a) Ising domain at array edge with weak bond at 
edge, (b) Edge domain for strong bond at edge. 

To check the above qualitative predictions, we have 
performed simulations on arrays of both geometries and 
with both uniform current injection and injection from 
superconducting busbars. Also, to avoid the effects of 
finite applied current which causes-a non linear IV rela- 
tion in the thermodynamic limitcrE3, we have computed 
the linear resistance from a linear response expression 
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in terms of equilibrium quantities. Our simulations are 
performed using the Langevin dynamical method of Falo 
et alB . We assume that the junctions on alternate 
columns have critical currents r]Io and Io elsewhere and 
are shunted by equal resistances R s . Each superconduct- 
ing grain has a small capacitance C to ground. The dy- 
namical equations for the phases Oi and the voltages V, 
of the grain at site i follow from charge conservation and 
the Josephson equation 

ddi/dt = 2eVi/h 

CdVi/dt = I E V13 sin(0j - Oi - Ay) 

<j> 

+R 7 ij2(v j -v i ) + j2nj (4) 

<j> <3> 

where rjij — 1 on all horizontal bonds and rjij — 77, 1 on 

alternating vertical bonds as in Fig .(1) and the sums over 

< j > are over the nearest neighbors of site < i > . The 

thermal noise current JH 1 on the bond < ij > is Gaussian 

i j ... 
distributed and obeys the fluctuation dissipation theorem 

< r$(t)I&tf) >= (2T/R s )(S ik S jl - S«6 jk )5(t - *') (5) 

For the columns of grains at either edge of the array 
labelled by i = 1, 2, • • • , L y 

CdVf/dt = 7 + 1 E Vij sm(0j - 0< - A tj ) 
<j> 

CdVf /dt = -J + J E Vij wx(6j -81- Mi) 
<j> 

+ j r; 1 E(^-^) + E / ^ ( 6 ) 

<i> <j> 

where I is the current per bond injected into each grain 
on the left column and extracted from each grain on the 
right column. For convenience, we choose system sizes 
with the number of columns L x an odd integer with free 
boundary conditions at the left and right hand edges and 
the number L y of horizontal rows an even integer with 
periodic boundary conditions in the vertical direction. 
We choose to use the Landau gauge where Afj = on all 
horizontal bonds, A?- = it on the odd numbered vertical 
columns (first, third, • ■ •, last) and A v - — on even num- 
bered columns. This is a particularly convenient choice of 
lattice size and gauge as it permits an integer number of 
2x2 unit cells and also for a simple formulation when su- 
perconducting busbars are connected to the left and right 
hand edges of the array by a set of L y junctions as the 
phase and voltage on each busbar is then y-independent. 
The equations governing these phases 0L(t), On(t) and 
the voltages V L (t), V R (t) are thenM 

dO L . R /dt = 2eV LM /h 

CdV L /dt = I + Ly 1 E [la sin(0; - L ) 
1=1 



+ R- 1 (V i -V L ) + I%] 

i—Ly 

CdV R /dt =-I + Ly 1 E [!o Mb - r) 

i=l 

+ R- l {V l -V R )+lt h R ] (7) 

Here the sums over i are over the L y sites connected to 
the busbars. The mean voltage drop V across the system 
is given by 

V/RJo = m r ) - <t>(P))/tr (8) 

where <p(t) — L (t) — 6 R (t) is the phase difference across 
the array at time t and t r is the run time of a simulation. 
The temperature T is measured in units of hIo/2e and 
the time t is in units of l/u>j = {TiC /2eIo) 1 ^ 2 , the inverse 
Josephson plasma frequency. A time step of typically 
At = 0.05 in these unitsca was used in the numerical 
integration. Changing the time step does not change 
the results. Most of our simulations have been done on 
small systems with L x = 17 and L y = 16 with periodic 
boundary conditions in the y direction. In any event, we 
found that boundary effects in the transverse y direction 
are not significant to the accuracy of our simulations. 

The IV relation is shown in Figs. (3) and (4) for the 
anisotropy parameter r\ = 0.5 at different temperatures 
T with the current I applied perpendicular to the weak 
bonds using the method of uniform injection described 
by Eq.(^|) with periodic boundary conditions in the trans- 
verse y direction. For the array of Fig. (2a), the onset 
of linear Ohmic behavior V = RlI is consistent with the 
data for T > 0.2 which is close to the estimate of the Ising 
critical temperature T c j rj 0.18 for this value of 77 = 0.5. 
When the geometry is changed so that the strong bonds 
are at the edges as in Fig. (2b), Ohmic behavior is ob- 
served only at higher temperatures T = 0.5—0.6, which is 
closer to the estimate of the XY transition at T c xy ~ 0.5 
as expected. For 77 — 2, as expected from the qualitative 
arguments, the situation is reversed and the results of the 
simulations are shown in Fig. (5) for arrays with the rj or 
strong bonds at the edges as in Fig. (2a) and in Fig. (6) 
for arrays of with the rj bonds not at the edges as in Fig. 
(2b). In the latter case, the onset of linear dissipation 
is observed at T « 0.4, close to T c j while in the former 
case at T ks 0.6, closer to the XY transition. Again this 
agrees with our qualitative arguments. 
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FIG. 3. IV characteristics for array with weak edge bonds 
as in Fig. (2a). r\ = 0.5, L x = 17, L y = 16 and periodic BC 
in y direction. 
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FIG. 4. IV characteristics for array with strong edge bonds 
as in Fig. (2b). 77 = 0.5, L x — 17, L y — 16 and periodic BC 
in y direction. 
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FIG. 5. IV characteristics for array with strong edge 
bonds. 77 = 2, L x = 17, L y — 16 and periodic BC in y di- 
rection. 
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FIG. 6. IV characteristics for array with weak edge bonds. 
77 = 2, L x — 17, L y — 16 and periodic BC in y direction. 

When the applied current / is finite, the IV relation 
always has a non-linear contribution. The resistance, 
defined by R — V/I, is proportional to I~ 7tK ^( t ) for 
T < T c j in the thermodynamic limit which may obscure 
the small predicted linear resistance of 0(1/L). To ob- 
tain a more definitive signature which is free of these 
non-linear finite current effects, one can define the linear 
resistance i?i by a Kubo formula in terms of an equilib- 
rium voltage- voltage correlation function at / = 

/+00 
dt < V(t)V(0) > (9) 
-00 

where V(t) is the total size-dependent voltage across the 
array. The results of the dynamical simulation using the 
Kubo expression of Eq.(||) at / = are shown in Fig. (7) 
together with those from a simulation with I/Iq = 0.02, 
the lowest accessible current, and defining Rl — V/I. 
This is always larger than the value obtained from Eq. @ 
because of the non-linear contribution when / > 0. This 
can be regarded as confirmation of our arguments that 
the onset of linear dissipation at T c i is not an artifact 
but a real effect caused by the Ising transition and may 
be observable by experiment. Our results for arrays of 
size L x = 17, L y — 16 with rj = 0.5 and periodic bound- 
ary conditions in the y direction are summarized in Fig. 
(8) where we show the linear resistance Rl(T) computed 
from the Kubo expression of Eq.(j^) for arrays geome- 
tries of Fig. (2a) and Fig. (2b). From this it is clear that 
there is an onset of linear resistance around T c j when 
the weak bonds are at the edges and at a higher tem- 
perature T f=a T c xy when strong bonds are at the edges. 
Our results for the phase diagram in the (T, rj) plane are 
shown in Fig. (9) where we plot T c j(rj) determined by 
purely equilibrium Monte Carlo simulations and also by 
the onset of Rl for arrays of Fig. (2a) from Eq.(^|). To 
the accuracy of our simulations, the two methods agree 
providing additional evidence to support our qualitative 
arguments. As a final check, we performed a series of sim- 
ulations at finite I on larger systems up to L = 200 to 
check the prediction R L — 1/L of Eq. (|| ) for T > T cI and 
in Fig. (10) we show VL/I for several values of L and we 
see that this is consistent with being an L-independent 
constant although the errors are fairly large. 
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FIG. 7. Linear resistance for array with weak edge bonds, 
r) = 0.5, L x — 17, L y — 16 and periodic BC in y direction, as 
estimated from the Kubo formula of Eq.^ at / = and from 
V/I at I = 0.02 




FIG. 8. Linear resistance Rl as a function of T from 
Kubo formula of Eq.(|^) at I — for array with 
L x — 17, L y — 16, rj = 0.5 in Fig. 2 with weak edge bonds 
(a) and with strong edge bonds (b). Note that Rl(T) > 
for T c xy > T > T cI for weak edge bonds and T > T c xy for 
strong edge bonds. 




FIG. 9. Onset temperature T for R L {T) > from Kubo 
formula Eq.(^) as function of rj denoted by solid dots. Ising 
and KT transition temperatures from equilibrium Monte 
Carlo simulations are open diamonds and pluses, respectively. 

The temperatures at which the onset of linear resis- 
tance occurs is consistent with the Ising critical temper- 
ature and happens for arrays of Fig. (2a) when r\ < 1 and 
Fig. (2b) when 77 > 1. This suggests that the linear resis- 
tance observed when T c j < T < T c xy is a result of both 
Ising disorder and free boundary conditions at the array 
edges allowing for motion of fractional charges induced 
by the domains, as predicted by qualitative arguments. 
The presence of weak bonds at the edges of the arrays 
of Fig. (2a) for 77 < 1 and of Fig. (2b) for 77 > 1 is not 
essential as changing the strength of the edge bonds does 
not change the scenario in any qualitative way but only 
quantitatively. The determining factor is that the peri- 
odicity of the Ising domain walls is two lattice spacings 
which is a consequence of the energy of a domain wall on 
a weak bond being larger than on a strong bond for rj < 1 
and the reverse for 77 > 1. When the edges correspond to 
a weak bond, domains at edges have width 71+ 1/2 unit 
cells which implies that there are free fractional charges 
at the ends of these domains in thermal equilibrium and 
these are free to move under the influence of an applied 
current. The strength of the edge bonds does not af- 
fect the argument in any essential way except to alter 
the magnitude of the dissipation. We have checked that 
the dissipation onset at T c j is in fact an edge effect by 
performing simulations in a busbar geometry where the 
current is injected and extracted by attaching the edges 
of the array to superconducting busbarsE£l described by 
Eq.(^). In this geometry, fractional charges (vortices) 
are repelled from the edges and domain formation is sup- 
pressed, thus minimizing or completely suppressing edge 
contributions to the dissipation. The IV relation for such 
an array of Fig. (2a) with 77 = 0.5 is shown in Fig. (11) 
for the current I normal to the weak bonds and in Fig. 
(12) for I parallel to the weak bonds. There is good 
agreement between the TV relations of Fig. (11) and 
Figs. (4,5) in which there are strong bonds at the edges 
and where we expect there to be no linear contribution. 
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One also observes that the slopes of the IV curves in the 
busbar geometry are independent of the direction of / as 
the slopes in Figs. (11) and (12) are the same for the same 
temperature T . This slope a(T) = log V/ log J is plotted 
in Fig. (13),-tpgether with that deduced from the helic- 
ity modulus^ 7 = (l x j y ) 1/2 byN a(T) = iry/T. The 
agreement between the different methods is by no means 
perfect but we believe the numerical support for our qual- 
itative predictions is more than adequate to demonstrate 
their validity. 
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FIG. 12. IV characteristics of array 

with L x = 17, L y = 16, r\ = 0.5 and current in y direction 
along weak bonds. 
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FIG. 10. Plot of (V/I)L at T = 0.4, / = 0.02 against array 
size I. 
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FIG. 11. IV characteristics of arrays with weak bonds at 
edges in busbar geometry of size I x — 17, I y = 16, rj — 0.5 
with current in x direction, normal to weak bonds. 
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FIG. 13. Exponent a(T) from IV curves and helicity mod- 
ulus. 

The conclusion we reach is that the low T Ising tran- 
sition in anisotropic arrays of Josephson junctions does 
have an effect which may be observable in zero frequency 
transport measurements. This effect is the onset of a 
linear resistance at T c j in arrays of the appropriate ge- 
ometry. Although this is an edge effect with Rl(T) ~ 
(T — T c j)/L it should be large enough to be measurable 
in an array of reasonable size. This onset is an unambigu- 
ous signature of the low T Ising transition as, in the ther- 
modynamic limit R L (T)L = A{T) ~ (T — T cI ) and bulk 
finite size corrections to this will decay very rapidly with 
increasing L and should be negligible compared to the 
X/L contribution discussed here. These bulk finite size 
effects can be estimated by performing the same mea- 
surement with the current parallel to the weak bonds 
when the 1/L contribution to Rl will be absent, which 
provides a method of assessing the feasibility of the pro- 
posed measurement. This small edge contribution to the 
linear resistance is a real effect which is a signal of the 
bulk Ising transition but, as is well known, the inter- 
pretation of such measurements is very difficult at tem- 
peratures well belojsaj'JJie KT temperature T c xy because 
of screening effectsoEj. We have not attempted to take 
such effects into account in this work but there is a pos- 
sibility that they may invalidate our conclusions. Clearly 
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this should be investigated as should various other effects 
such as disorder in the junction strengths which will in- 
evitably be present in a real system, vortex and domain 
wall pinning etc. It is not very aesthetically satisfying 
that one must resort to an edge effect detectable by a 
voltage measurement at finite current as there are limits 
to the sensitivity_of this, especially because a typical ex- 
perimental arrayO is much larger than those of this work 
so that the dissipation at the edges will be much smaller 
than in our simulations. In fact, in the thermodynamic 
limit L — > oo the linear dissipation in the temperature 
range T c j < T < T c xy where the vortex lattice is melted 
should vanish. It would be much more satisfying if some 
equilibrium bulk quantity would give a signal of the Ising 
transition butriHpview of the very weak signal in the he- 
licity modulual3'ta, this will be a very difficult effect to 
observe. There is the possibility that .a flux noise mea- 
surement similar to that of Shaw et alB may show a de- 
tectable signal of the Ising transition. However, naively 
this seems not very hopeful as the noise spectrum basi- 
cally measures the time dependent correlation function 
< N(t)N(0) > where N(t) is the total vorticity (charge) 
enclosed by the squid detector. Since the Ising transition 
is signalled by the proliferation of domains with frac- 
tional corner charges but with zero net charge, as in Fig. 
(lb), it is difficult to see how these can contribute signif- 
icantly to the correlation function determining the noise 
spectrum. However, this is still worth investigating as 
it is one of the few remaining possibilities to detect this 
elusive transition. 
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